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We analyze the influence of the evolution of light absorption by gray dust in the host galaxies 
of type la supernovae (SNe la) and the evolution of the mean combined mass of close-binary 
carbon-oxygen white dwarfs merging due to gravitational waves (SNe la precursors) on the 
interpretation of Hubble diagrams for SNe la. A significant increase in the mean SNe la energy 
due to the higher combined masses of merging dwarfs should be observable at redshifts z > 2. 
The observed relation between the distance moduli and redshifts of SNe la can be interpreted 
not only as evidence for accelerated expansion of the Universe, but also as indicating time 
variations of the gray-dust absorption of light from these supernovae in various types of host 
galaxies, observational selection effects, and the decrease in mean combined masses of merging 
degenerate dwarfs. 
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1 INTRODUCTION 

This paper continues our study [Ij, where we demonstrated that the mean combined 
masses of binary carbon-oxygen white dwarfs (CO WDs) merging due to gravitational- 
wave radiation should evolve with time. Paper [T] presents the results of population 
syntheses performed with the "Scenario Machine" (a computer code for evolutionary 
studies of close binaries), as well as computations of the merger rate and mass distributions 
of merging white dwarfs (WDs). Such studies are important because observations of 
distant type la supernovae (SNe la) in the late 1990's led to the conclusion that the 
expansion of the Universe was accelerating^ ^ Ej. Thus, understanding the physics of 
SNe la determining the possible evolution of their brightness with time is needed to 
understand the fundamental properties of the Universe. In addition, in contrast to our 
first study, we also consider here the influence of the evolution of gray dust in the host 
galaxies of SNe la on changes in the absorption of light from these supernovae. Light 
absorption by gray dust is able to produce the same effect of fainter distant supernovae 
as the proposed accelerated expansion of the Universe. 



^At z « 1 were found to be fainter than expected for standard supernovae of this type by, on average, 
about 0™.2 [113]. 
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There are two main models for the origin of a SNe la outburst. The first requires a close 
binary consisting of a CO WD and a non-degenerate star, whose matter is accumulated 
onto the white dwarf via accretion and hydrogen and helium burning in shell sources 
in the dwarf's envelope. As soon as the WD mass exceeds the Chandrasekhar limit, a 
thermonuclear explosion occurs (or an explosion of the helium shell causes detonation 
of carbon in a dwarf with a mass below the Chandrasekhar limit), observable as a SN 
la The second model requires a close binary consisting of two CO dwarfs. Under 
the influence of gravitational- wave radiation, the binary's semi-major axis decreases, the 
dwarfs merge, and the resulting thermonuclear explosion is manifest as a SN la [SI E]^- 
We consider the second scenario to be the main one 

Yungel'son [7] aims to prove that mergers are the principal, and possibly the only, 
source of SNe la. The accretion scenario is probably unable to explain SNe la: CO WDs 
are probably not able to accumulate enough gas to make their masses exceed the Chan- 
drasekhar limit [7J. Observations also present evidence in favor of the WD-merger scenario 
[H E]. However, the viewpoint of [3 |9] are not shared by all authors [10]. Nevertheless, 
the main mechanism for SN la explosions in elliptical galaxies with ages in excess of ap- 
proximately 200 million years is CO WD mergers [9]- [12]. The energies of some SN la 
correspond to masses of burned carbon and oxygen 2Mq [T3], providing direct evidence 
in favor of the dwarf-merger model and confirming our estimates [T]. 

The merger rate of close binary WDs and its dependence on the star-formation history 
were studied in [Ullll]; the results of later studies are not very different (see, for example, 
[T^ [15]). This may be because the WD merger rate is virtually independent of the 
merger mechanism if the merger time scale depends only on the initial semi-major-axis 
distribution for binary WDs [TU]. Merger candidates are actively being searched for, and 
merger rates and the delay of the merger of degenerate close-binary components estimated 
(see, for instance, [T6| fT7]). 

The role of the possible evolution of the energy of SNe la was studied ten years ago in 
|18j . where several possible hypothetical luminosity evolutions for SNe la were considered. 
According to those computations, introducing even a small amount of evolution makes it 
very difficult to determine cosmological parameters such as the dark-energy density. We 
demonstrated in P that SNe la can be considered standard candles only after approxi- 
mately one billion years of galaxy evolution^. In the course of their evolution, the mean 
energy of SNe la should decrease by approximately 10%, and the difference between the 
highest and lowest energies of SNe la should be at least a factor of 1.5"^ at both high and 
low redshifts. These obvious observational selection effects must be taken into account 

^In this case, the combined mass of the dwarfs probably exceeds the Chandrasekhar hmit. 

■'This result is independent of the star-formation history for the three versions of the star-forming func- 
tion considered in [1] , all stars are born simultaneously (star formation in an idealized elliptical galaxy) , a 
constant star- formation rate (star formation in an idealized spiral galaxy), and the star formation history 

of ng. 

^It is a factor of two in the extreme case: a merger of two WDs with a combined mass approximately 
at the Chandrasekhar limit, or a merger of two WDs, each having a mass close to the Chandrasekhar 
limit. However, mergers of CO WDs with combined masses in excess of 2.1Mq are rare. In addition, 
if we take into account the hypothetical possibility of explosions from CO WD mergers with combined 
masses below the Chandrasekhar limit, the difference between the highest and lowest energies of SNe la 
could be even greater than a factor of two. Studies of possible sub-Chandrasekhar scenarios for the origin 
of SNe la have now begun (see, for instance, [501 [H]). 
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when estimating the parameters of the Universe's expansion using SNe la. The maximum 
brightness of such supernovae is not a universal constant, and has some dispersion. While 
we see "all" the closest SNe la, only the brightest are detected at increasing distances. 
This could imitate a brightening of SNe la with distance. The magnitude of this effect is 
comparable to the intrinsic brightness dispersion for this type of supernovae. 

The brightness of a SNe la can be attenuated by absorption in matter around the 
supernova precursor, as well as absorption by dust in the supernovas parent galaxy; the 
principal characteristics of dust in the host galaxies of supernovae, total dust mass, dis- 
tribution of dust in the galaxy, and dust composition can all vary in time. Redshift- 
dependent absorption by intragalactic gray dust could introduce errors into the extent 
to which SNe la can be considered standard candles; not taking into account absorption 
by gray dust makes photometric distances too large compared to a dust-free Universe 
and could give rise to an apparent "acceleration" of the expansion of the Universe [22]. 
According to ^23j , taking into account absorption by gray dust can make a SNe la fainter 
by O.^OS at a redshift of z = 1.7, potentially capable of causing disagreement between 
Hubble diagrams calculated assuming a dark-energy density fl\ ~ 0.7 and observations 



Absorption by intergalactic gray dust was suggested in [211 [25] as an explanation of 
the fainter magnitudes of supernovae at large distances, as an alternative to accelerating 
expansion of the Universe, almost immediately after the classic papers [21 13] on observa- 
tions of distant supernovae. In addition, a model with intergalactic dust was considered in 
[26] . where it was suggested that intragalactic gas was not gray, and so could be accounted 
for by means of precision spectroscopy. However, later studies gave somewhat different 
results 1271. 



A comparison of the brightnesses of some 85 thousand quasars and the positions of 
about 20 million galaxies (based on the SDSS survey) demonstrated the presence of dust in 
intergalactic space, comparable in amount and composition to the dust in galactic disks 
[28] . Problems related to the correction of supernova parameters taking into account 
absorption by this dust are discussed in [21], where the influence of gray dust is not 
considered. The authors of [2Sl 121] note that if the detected dust does not contain gray 
dust, they are able to provide an upper limit for absorption, but if gray dust is present, 
their computations give a lower limit for the influence of dust. 

The example of the Milky Way shows that absorption in the parent galaxy can reach 
1"^, even in the polar direction. It also makes supernovae in disk galaxies fainter. For 
example, observations of stars in the Sloan Digital Sky Survey (SDSS) led to the discover 
and investigation of three absorbing clouds in the Galaxy, each hundreds of parsecs in 
size. The absorption in these clouds is 0™.2 — 0"\4, and the clouds exhibit gray, almost 
wavelength-independent absorption in the spectral range of the survey [27]. When esti- 
mating the evolution of the brightness of SNe la in time, we must remember that the 
gas component of the galaxies containing these supernovae also evolves. This introduces 
additional uncertainty into brightness estimates for supernovae. The higher infrared lu- 
minosities of distant {z 1) galaxies indicate the presence of high dust contents in young 
galaxies (e. g., [30]). 

We can put the following limit on radiation absorption by dust in an elliptical galaxy 
with a mass comparable to that of the Galaxy. Let the total mass lost by red supergiants 



of SNe la. 
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be ~ IMq per year. The mass loss time scale for an elliptical galaxy is ~ 10^ years. Thus, 
the stationary amount of gas in the galaxy is ~ lO^M©. In this case, the absorption for 
the usual abundance of dust will be ~ 0™".!, if we assume the same efficiency for dust 
absorption in elliptical and disk galaxies. 

To estimate the absorption of light from supernovae in disk galaxies, we applied the 
evolutionary model of [31]- [33]. The assumed mass and radius of the model galaxy were 
2 ■ IO^^Mq and ~ 10^ pc. We estimate the star-formation rate from the condition that the 
gaseous disk be fully ionized by young massive stars [31]. The thickness of the gas disk is 
estimated from the condition that the rate at which energy is supplied to the turbulent 
motion of the galaxys gas by supernovae (II and Ib,c) be the same as the dissipation rate 
of the turbulent motion of this gas |3lj. In this case, the rate of supernova outbursts 
is determined by the star-formation rate in the model galaxy, and the dissipation of the 
turbulent motion of the interstellar gas is determined by collisions between interstellar 
gas clouds. 

Modeling demonstrated that the star-formation rate in a disk galaxy varies consider- 
ably in time [32]. The star-formation rate increases during the several first billion years. 
In a galaxy similar to ours, the star-formation rate reaches a maximum of ~ IO^Mq per 
year at an age of 2 — 4 billion years. During the subsequent evolution, the star-formation 
rate gradually decreases to its current level. The model takes into account gas enrichment 
with heavy elements, making it possible to follow the time evolution of the gas compo- 
nents optical depth when checking the mass of the gaseous component. Knowing this, 
we are able to estimate the time evolution of the apparent brightnesses of supernovae, 
assuming that all are located in the planes of their host galaxies, which are orthogonal 
to the line of sight. In this way, we can take into account the effect of absorption of 
light from supernovae in disk galaxies. For the current epoch, this amounts to ~ O^.S, 
and increases for young galaxies with high gas and dust contents. Dust is thus able to 
attenuate the brightness of SNe la. 

It is especially important to note that we are concerned here with gray dust, whose 
absorption is wavelength- independent. The existence of such dust in the galactic plane 
is now beyond doubt [27j. Fine dust causes reddening of light and thus, in principle, 
can be taken into account. It should also be borne in mind that modeling of the time 
evolution of a galaxy's optical depth remains uncertain, and the optical depth itself cannot 
be estimated better than to within a factor of two [311 [32] . Gray-dust light absorption 
occurs independently of the supernova model that is preferred. We will assume here that 
the amount of gray dust in a parent galaxy of a SN la is proportional to its total amount 
of dust. 

2 MODEL EVOLUTION OF SNe la BRIGHTNESS 
WITH INCREASING REDSHIFT 

In this study, we analyze the computational results used to plot Figs. 1-3 in [Ij^. We 
used these computations to derive the mass distributions of merging WDs [1] (Fig. 1) and 

^The computations of [T] were based on the "Scenario Machine". A description of the code can be 
found in [341 135) . and of the population synthesis, in |36| . 
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the dependence of the mean mass of merging WDs on the time elapsed since the onset 
of star formation [Ij ( Fig. 2). Virtually independent of the form of the star- formation 
functions, the mean mass of merging CO WDs with combined masses in excess of the 
Chandrasekhar limit decreases from 1.9Mq to ~ 1.7Mq within about one billion years 
after the onset of star formation, after which the mean mass of merging WDs remains 
practically the same. 

The aim of this study is to provide a theoretical analog of a Hubble diagram p] (Figs. 
4, 5) taking into account the brightness evolution of SN la and the evolution of dust 
absorption, using the theoretical model presented in [33j. We assume that the absolute 
magnitude of a SN la depends on the mean mass of the merging CO WDs as 

M = C-2.5-logMs, (1) 

where M is the absolute magnitude of the SN la, the combined mass of the merging 
CO WDs, and C is a constant. Thus, the higher the combined mass of the merging 
dwarfs, the brighter the supernova. Using [1] to re-calculate the data plotted in Fig. 2 of 
[Ij, we obtained the dependence of the absolute magnitude of a SN la on the time since 
the onset of the formation of the corresponding stellar population. 

The resulting dependence of the absolute magnitude on the time from the onset of the 
galaxy's evolution until the supernova outburst must be transformed into the dependence 
of the observed SNe la magnitudes on the observed redshifts for the host galaxies. These 
computations used Eq. 20 (to calculate the time-redshift relation) and 28 (to calculate the 
photometric distance- redshift relation) from [37]. The free parameters in these formulas 
are the density of dark energy Qx, the matter density Qm, the total density of the Universe 
^Itot, and the current Hubble constant Hq. Another free parameter of our problem is the 
redshift z* corresponding to the onset of star formation in galaxies. 

The last stage of our computations is to calculate the dependence of A(m — M) on 
redshift z: 

A(m - M) = (mi - C) - (ma - C) = mi - m2, (2) 

where m and M are the apparent and absolute magnitudes, mi the apparent magnitude 
calculated taking into account the time evolution of the mean mass of the merging WDs, 
and m2 the apparent magnitude calculalated assuming that a SN la is a standard candle 
with absolute magnitude C. 

We took the dependence for light absorption by gray dust from [33j (Fig. 2 presents 
the dependence of on the redshift z; we used the curve corresponding to a disk galaxy 
with the mass-radius relation M ~ R^, see [33] for details). The redshift dependence of 
the absorption in this figure was derived using the formula t{z) = 3 jj-^^ (1^^)3/2 relating the 
time t{z) elapsed since the Big Bang and the redshift z. We place this relation on the 
Hubble diagram in order to qualitatively illustrate the evolution of dust absorption in 
disk galaxies, compared to the accelerated fading of supernovae at redshifts 0.5 — 1. Note 
also that absorption by dust grains of any size was considered in [3,3] ; that paper contains 
no estimate of the fraction of gray absorption due to grains with sizes in excess of l/im. 
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3 RESULTS 



The results of our computations are presented in Figs. [Hand [21 The values computed 
using [2] are plotted versus redshift. The magnitude m2 used in [2] to calculate the curves in 
Figs. [Hand [2] is the dependence of the apparent magnitude m on the redshift 2, assuming 
that all SN la are standard candles with absolute magnitude C and that they do not 
undergo any absorption (at least no time- variable absorption). This relation corresponds 
to the horizontal straight line at zero magnitude (curve 3). The magnitude mi is the 
dependence of m on ^ taking into account one of the studied effects. 

Curve 1 describes the dependence of A(m — M) on z, assuming that either absorption 
plays no role or it is constant in time, and that the absolute magnitude evolves in time 
from its maximum value at the beginning of the evolution to its minimum, approximately 
constant, value after about a billion years of stellar evolution in galaxies. The adopted 
cosmological parameters are VL\ = 0, = 0.2, Qtot = 0.2, Hq = 70 km s~^ mpc~^, 
z* = 10. 

Curve 2 describes the dependence of A(m — M) on z, assuming that evolutionary 
effects related to the mass of the merging dwarfs and light absorption in galaxies play no 
role; here, we adopted Qx = 0.7, Qm = 0.3, Qtot = 1, -f^o = 70 km s~^ mpc~^. 

The line 3 corresponds to a model without evolutionary effects and with no influence 
from dark energy {Qx = 0, Qm = 0.2, fltot = 0.2, = 70 km s~^ mpc~^). 

Curve 4 describes the evolution of light absorption by dust. This relation was taken 
from [33] (Fig. 2 displaying the dependence of the absorption on the redshift z, 
adopting the curve corresponding to a disk galaxy with M ~ i?^). 

Figure [T] shows that dust absorption is able to attenuate the brightness of SN la by as 
much as ~ 1™.3. This absorption is reached at redshifts z ~ 3. The evolution of the mean 
mass of the merging CO WDs begins to be appreciable at ^ > 2 and leads to an increased 
brightness of SNe la in the past compared to the current epoch. The mean apparent 
magnitude can decrease by ~ 0^.2 at 2 > 8, in full agreement with the possibility that 
the mean mass of merging CO WDs (and, according to our assumptions, the mean energy 
of SNe la) was higher at the beginning of galaxy evolution, in the distant past, compared 
to the current epoch. 

To make it easier to compare our results to Figs. 4 and 5 from [^ , Fig. [2] isolates 
the redshift range between 0.01 and 1. The evolution of the mean mass of the merging 
dwarfs essentially cannot influence the results in this redshift range. Dust absorption 
attenuates the brightness of supernovae much more strongly than if we allow for a possible 
contribution from dark energy (Figs. 4 and 5 in [2]), but recall that [33] presents data for 
absorption by any dust, including dust whose influence can be taken into account based 
on reddening. In order to explain the accelerated fading of SNe la at ^ = 0.5 — 1 as due to 
evolution of gray-dust absorption in disk galaxies, it is sufficient to suppose the amount 
of gray dust is proportional to the total amount of gas in the host galaxies of supernovae, 
and comprises several tens of percent of the total dust absorption in these same galaxies. 
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4 CONCLUSIONS 



We have considered two possible channels for the "standard" behavior of the SNe la 
brightness being violated. The first is the time evolution of the mean mass of merging CO 
WDs, which are the precursors of SNe la. The second is the time evolution of absorption 
by gray dust in the SNe la host galaxies. 

According to our computations, the mean energy of SNe la should increase from z > 2 
and become significantly higher at 2; > 8 due to the fact that, on average, mergers involved 
higher-mass WDs during the early evolutionary stages of the Universe than at the current 
epoch. Such distant supernovae have not yet been discovered. Thus, the conclusion that 
the expansion of the Universe is accelerating based on observations of SNe la at redshifts 
to 2; ~ 1 made in [21 E] will not be affected by evolution of the mean mass of the merging 
dwarfs. At the same time, it should always be borne in mind during such estimation 
that our models, which are based on population syntheses, assume that the objects whose 
evolution we are studying (galaxies) were born at some time in the past (in the current 
study, at z* = 10), then evolved. However, if a burst of star formation occurred at some 
epoch after z*, conclusions about the increase of the mean SNe la energy should be related 
to the time since that burst. Independent of the population's age, there will exist some 
dispersion in supernova parameters, because the combined mass of the merging dwarfs 
can vary within a factor of approximately 1.5. 

According to the model considered here, the evolution of radiation absorption by 
galactic dust can give rise to a much stronger attenuation of the brightness of distant 
supernovae than the assumption that dark energy is present in the Universe. At the same 
time, the absorption described by curve 4 in Figs. 1 and 2 is the sum of the absorption 
that can, in principle, be taken into account based on reddening and gray absorption, 
which cannot. By varying the model using its intrinsic parameters and introducing the 
required fraction of gray dust, we are able, in principle, to achieve an essentially ideal 
agreement between curves 2 and 4 over a wide redshift range (this reasoning is very similar 
to the considerations presented in [21]- [2S])- However, this approach is often criticized: 
fine tuning of the model parameters is believed to be able to reproduce any results, and 
interpretation of the phenomenon in this way seems unattractive [38] . 

During the analysis of the SDSS-II supernova survey (and of the light curves of these 
supernovae), SNe la were found to be approximately 0^.1 brighter in quiet galaxies than 
in galaxies with active star formation [39 f. The characteristics of absorption by dust are 
also different in galaxies with and without active star formation [39]. A dependence of 
the SNe la luminosity on galaxy type was also noted in [iO] . 

We suggest the following interpretation of our results. Figure 2 shows that absorption 
by dust in disk galaxies grows at essentially the same rate as the rate of dimming of 
supernovae due to accelerated expansion of the Universe. Hence, we suggest that the 
decrease in the observed brightness of supernova (on average, and this is important) is 
due primarily to radiation absorption by gray dust, whose influence cannot at this time 
be taken into account using standard procedures for correcting for absorption. However, 
further increase of the redshift results in increasing absorption, causing still stronger 
brightness attenuation of supernovae compared to a model without dust. At the same 

^Based on analysis of supernovae with redshifts up to z ~ 0.2 - see |39], Fig. 3. 
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time, some distant supernovae are found on the curve corresponding to decelerating, rather 
than accelerating, expansion of the Universe (see, for instance. Fig. 10 in [41j). This could 
be considered unequivocal evidence against the idea that the observed supernova fading 
can be explained by time-evolving absorption by gray dust^ if all host galaxies of SNe la 
were disk galaxies with active star formation. However, this is not the case: for example, 
the parent galaxy of SN 1997ff (one of the most distant supernovae, used as a standard) 
is elliptical [12] . 

As noted above, compared to disk galaxies, absorption in elliptical galaxies is insignif- 
icant, and apparently does not evolve over a considerable time interval, comparable to the 
Hubble time. Note also that population syntheses of the evolution of close binaries (e.g.. 
Fig. 1 in [12] or Fig. 2 in [15]) show that the rate of SNe la outbursts in an elliptical 
galaxy decreases by approximately a factor of two between 1 and 10 billion years after 
the formation of the elliptical galaxy. This means that the probability of observing a SNe 
la in an elliptical galaxy, rather than a disk galaxy, was much higher in the past (among 
supernovae with the highest z) than at the current epoch (for low z). This can explain 
the observed lack of "accelerated" supernovae for redshifts z > 1. 

Thus, when using SNe la as cosmological standard candles, it is necessary to take into 
account not only model characteristics of the supernovae, but also evolution of absorption 
by dust, including gray dust, in the supernova host galaxies. It is preferable to use super- 
novae in galaxies without star formation as distance indicators^. Supernova outbursts in 
galaxies with active current star formation can serve as probes for studies of the evolution 
of absorption in such galaxies (including gray absorption^). 

Variations in the combined masses of binary degenerate dwarfs with ages comparable 
to the age of the Universe do not significantly change the brightness of SNe la explosions. 
It follows from the discussions above that analysis of the evolution of SNe la brightnesses 
aimed at improving cosmological models must be carried out separately for elliptical and 
spiral galaxies. SNe la in elliptical galaxies seem to provide a more reliable cosmological 
tool than SNe la in spiral galaxies, since the latter contain dust. Searches for and studies 
of very distant SNe la, at redshifts z > 2 are especially important. 
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FIGURE CAPTIONS 

Fig. [H Difference between the magnitude-redshift dependences taking into account 
the effects studied here (time evolution of the absorption and the mean mass of merging 
WDs) and assuming standard SNe la that are not subject to variable absorption in time, 
A(m — M). Curve 1 assumes that absorption does not play any role and that the absolute 
magnitude evolves in time from a maximum value at the beginning of the evolution to a 
minimum, approximately constant, value after about one billion years of stellar evolution 
in galaxies (for = 0, = 0.2, = 0.2, Hq = 70 km s"^ Mpc-\ z* = 10). Curve 2 
assumes that absorption does not play any role and the absolute magnitudes experience 
no evolutionary effect {Qx = 0.7, Qm = 0.3, Qtot = ^, Hq = 70 km s~^ Mpc~^). The line 
3 corresponds to a model without evolutionary effects {Qx = 0, Q^j = 0.2, Qtot = 0.2, 
Hq = 70 km s~^ Mpc~^). Curve 4 represents the evolution of light absorption by dust, 
taken from [3^ (Fig. 2 displaying the dependence of on the redshift z; the curve 
corresponding to a disk galaxy with M ~ i?^). 

Fig. 121 Same as Fig. 1 for redshifts between 0.01 and 1. 
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